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We analyze a quantum mechanical gyroscope which is modeled as a large spin and used as a 
reference against which to measure the angular momenta of spin- 1/2 particles. These measurements 
induce a back-action on the reference which is the central focus of our study. We begin by deriving 
explicit expressions for the quantum channel representing the back-action. Then, we analyze the 
dynamics incurred by the reference when it is used to sequentially measure particles drawn from a 
fixed ensemble. We prove that the reference thermalizes with the measured particles and find that 
generically, the thermal state is reached in time which scales linearly with the size of the reference. 
This contrasts a recent conclusion of Bartlett et al. that this takes a quadratic amount of time 
when the particles are completely unpolarized. We now understand their result in terms of a simple 
physical principle based on symmetries and conservation laws. Finally, we initiate the study of the 
non-equilibrium dynamics of the reference. Here we find that a reference in a coherent state will 
essentially remain in one when measuring polarized particles, while rotating itself to ultimately align 
with the polarization of the particles. 

PACS numbers: 



I. INTRODUCTION 



Physical systems are described in terms of kinematical 
observables such as position, momentum, and spin. In 
the presence of symmetries, only a subset of those quan- 
tities - the physical observables - are physically relevant, 
as they too must obey the same symmetries. One way of 
obtaining physical observables is to use some systems as 
references relative to which the properties of other sys- 
tems can be measured, yielding observables which are 
relational. The state of the reference on its own cannot 
be measured: only its correlations with other systems 
acquire physical meaning through direct measurements. 
In many circumstances, the existence of a suitable ref- 
erence is conspicuous, often blurring the distinction be- 
tween kinematical and relational observables. 

For example, the walls of a laboratory can be used to 
align a Stern-Gerlach apparatus to very high precision. 
This reference enables one to measure any component of 
an electron's spin, despite the rotational invariance of the 
underlying theory. While this statement may seem con- 
tradictory, it is only due to the fact that the reference 
was omitted from the description of the measured quan- 
tity. Formally, the spin of the electron is only determined 
relative to the orientation of the measurement apparatus 
(and ultimately the laboratory walls, the stars, etc.) and 
all symmetries are respected. When the reference sys- 
tems are classical, whether or not to include them in a 
given physical description is purely an aesthe tic choic e 
which has no impact on physical predictions |BRS05a |. 
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These "two sides of the same coin" are analogous to doing 
calculations in electromagnetism with a special choice of 
gauge or using the gauge invariant language: each ap- 
proach has its advantages. 

However, if the measuring device (or at least certain 
essential parts of it) are inherently quantum mechanical, 
these distinctions become more relevan t. For ex ample, 
it was pointed out long ago by Wigner Wig52a( that a 
quantum mechanical conservation law places constraints 
on the ability to obtain precise measur ements of certain 
observables (see also the strengthening [AY60a(). Similar 
arguments can be made for clocks, rods and gyroscopes, 
as they are also often taken for granted and omitted from 
the theoretical description of the situation at hand. 

In this paper, we study a model of quantum measure- 
ment in which the task is to measure the angular mo- 
mentum of spin-i particles along some given direction, 
in the presence of a rotational symmetry. The informa- 
tion about the axis along which to measure is supplied 
by a quantum mechanical reference system. We will not 
only study the limitations this model places on perform- 
ing perfect measurements - rather, we will focus on how 
the measurements affect the reference, and thus how they 
change the character of measurements performed with 
the same reference at later times. 

While omitting the description of the reference at hand 
is often entirely justified from a practical standpoint, this 
issue has been the source of a great deal of confusion and 
heated debates. One illustration, particularly disputed 
in the Bose-Einstein Condensation (BEC) community, 
is the definition of an order parameter as either a local 
expectation or a correlation function. The former idea 
leads to the concept of spontaneous symmetry breaking 
[And72a| . where the field acquires a finite expectation 
value prohibited by the symmetry. The latter idea views 
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one part of the BEC as a reference relative to which the 
expectation of the field at another point in the BEC can 
be measured, t hus illus trating; the importance of correla- 
tion functions [LegQlal ] . Similar debates have ta ken place 
with regard to the phase of a laser field [Mol97al | being ei- 
ther absolute or only defined relative to a reference pulse. 

Another important example occurs in general relativ- 
ity, where the group of symmetries consists of diffco- 
morphisms of spacetime. Indeed, the issues concerning 
physical observables in general relativity are among the 
most profound in physics as they raise, among other 
things, question s regarding the measurability of time 
AB61a, PW83b], as well as that of the spacetime metric 



Ber61aj itself. Once again, some of these conceptual dif- 



ficulties can be circumvented by a careful analysis of the 
role of t he reference systems invo lved in our description of 



nature (PW83bLlAK84allRov91a| . We expect that resolv- 
ing such foundational issues would be an important step 
towards the much-sought-after marriage between quan- 
tum mechanics and general relativity at the heart of cur- 
rent theoretical physics research. 

Quantum mechanical descriptions of reference systems 
have revealed important physical consequences. These 
includ e the possibility o f breaking certain superselection 
rules |AS67al . IKMP04 a] and the existence of an intrin- 
sic decoherence in quantum gravity |Mil91al . lGPP04al 
lPou06aL IMP06 aj. However, it is the study of quantum 
reference systems as physical resources that has recently 
received the most attention. Indeed, using the tools of 
quantum information science, one can quantify the cost 
of esta blishing a reference between two distant parties 
[PS01a |. the cost of realizing a reference with a d esired 
degree of accuracy [BRS04aj . etc. (see [BRS06a] for a 
recent review). 

The latest addition to the theory of quantum reference 
systems is the study of their degradat ion under repeated 
measurements |BRS T06al . lBRST06b| . Indeed, the very 
act of using a reference to measure a system typically 
decreases the accuracy with which it can make a subse- 
quent measurement on another system. Stud ying specific 
models, the authors of BRST06l lBRST06bj ' came to the 
conclusion that this measurement accuracy typically de- 
creases as the square-root \ft of the number of measured 
systems t, or in other words, that the so-called longevity 
of a reference is proportional to the square of its size. 
Here, we revisit one of their models and analyze its dy- 
namics in more detail. We show how their result follows 
from simple physical principles based on symmetries and 
conservation laws which are applicable even to general 
systems. We also find that the case analyzed by those 
authors is in fact quite singular and does not reflect the 
general behavior of a reference's longevity, which we in- 
stead find to scale linearly with its size in the generic 
case. As we observe in Sec. I IV A[ this distinction ex- 
ists because the source particles cause the reference to 
undergo a random walk whose dynamics are dominated 
by fluctuations when the source particles are unpolarized 
and by drift when they are polarized. 



The importance of studying quantum reference frame 
degradation stems not only from its repercussions on 
fundamental physics outlined above, but also from its 
relevance to realistic experimental settings (see e.g. 
[WDL+04aLlRBMC04 aj for a description of experimental 
realizations relatively similar to our theoretical idealiza- 
tion). In particular, such effects are likely to play an 
important role in the realization of quantum information 
processors. We will elaborate on this aspect in Sec. IVTl 

The rest of this paper is organized as follows. Sec. |TT] 
sets the stage by introducing the model that we study and 
describing the task it is meant to accomplish. Sec. IIIII 
presents an analysis of the asymptotic behavior of the 
reference system. The following section studies the ther- 
malization time required to reach this asymptotic state. 
In Sec. [Vj we derive semi-classical equations of motion 
that model the dynamics of the reference system during 
the equilibration process. The last section summarizes 
our results and discuss some of their implications. 



II. THE MODEL 

Our primary object of study is essentially a gyro- 
scope, i.e. a physical system that singles out a particu- 
lar direction in three-dimensional space. In principle, an 
ideal (classical) gyroscope would allow an experimenter 
to measure the spin of a spin-i source particle along the 
axis of rotation of the gyroscope. The quantum analog 
of a gyroscope is a system with a large amount of spin; 
we henceforth refer to such a system as a reference. The 
state of the reference is described by a density matrix p 
on the spin-£ irreducible representation (irrep) of SU(2). 
Throughout, we will consider £, or equivalently the di- 
mension d — 2£ + I, as an indication of the size of the 
reference. 

Ideally, the reference will be in a coherent state and 
will thus have a maximal amount of angular momentum 
concentrated in some direction n. We only allow the ex- 
perimenter to operate on the system and the reference 
in a rotationally invariant manner. Thus, the only non- 
trivial observable to which the experimenter has access is 
the total angular momentum of the source particle and 
reference. When used in this fashion, the quantum ref- 
erence differs from its classical counterpart in two funda- 
mental ways: 

1. The measurements will only be an approximation of 
what would be obtained by using the corresponding 
classical reference. 

2. Each time the reference is used, it suffers an in- 
evitable "back-action" which ultimately changes 
the character of future measurements. 

In any practical scenario, this reference would in gen- 
eral consist of a composite system, e.g. composed of the 
many electron spins of a ferromagnet. That only the 
SU(2) degrees of freedom represented by p are of inter- 
est implies that we are assuming that all the constituent 
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particles of the gyroscope are used in a symmetric man- 
ner, and so can be considered as a single particle with 
large angular spin £, such as a Rydberg atom. Other in- 
teractions will typically lead to a shorter longevity of the 
reference, but we are here interested in the fundamental 
upper bound which physics places on its longevity. 

While the physics we analyze is independent of any 
background reference frame, our mathematical analysis 
requires us to choose one. This artificial background co- 
ordinate frame should be thought of as a choice of gauge 
with respect to which the underlying physics is insen- 
sitive. In terms of these background coordinates, the 
SU(2) generators of the reference are kinematic observ- 
ables which we denote L = (L x , L y , L z ). These operators 
satisfy (see e.g. [Sak94al |) 



[Li, Lj] = ie i3k L k 



(1) 



where e is the completely antisymmetric tensor for which 
1 . The corresponding raising and lowering opera- 



xyz 



tors L± = L x ± iL y satisfy 
[L+,L_] = 2L Z , 



[L Z ,L ± ]=±L 



±- 



(2) 



Consider now a source which emits spin-^ particles, 
each described by the same density matrix £. The corre- 
sponding angular momentum operators S = (S x , S v , S z ) 
also satisfy the commutation relations given by Eq. |T]). 
The constrained relational observable allowing us to mea- 
sure the spin of the particle relative to the quantum refer- 
ence is the square of the total spin operator J 2 = (L+S) 2 . 
Since L 2 and S 2 are in the center of the algebra, J 2 (and 
arbitrary functions thereof) is the only non-trivial oper- 
ator that can be measured in our toy universe. Mea- 
surement of this observable projects the combined state 
of the system and the reference onto an eigenspace of J 2 
with eigenvalue + where j = £± |-. After this mea- 
surement, we assume that the measured source particle is 
discarded and that the measurement result is forgotten; 
in other words, the measurement scheme is non-adaptive. 
While it is an interesting question whether or not the 
history of results from subsequent measurements could 
potentially be used to give a more accurate description 
of the dynamics of the reference, we do not consider such 
approaches in this paper. 



A. The approximate measurement 

We now describe the measurement induced on the 
source particle by the joint measurement of J 2 . This 
projective measurement on the combined system induces 
a positive operator valued measurement (POVM) on the 
source. To find an expression for this POVM, we begin 
by writing the respective projections onto the j = I ± ^ 
irreps as 



± 2 I d 



(3) 



The measurement {A^ } induced on the source is now 
given by a partial trace over the reference, which yields 



A 



Tr R Il + {p®I 2 



h± 



4(L) ■ S + 1, 



Note that the measurement depends only on the expec- 
tations (Li) = Tr pLi of the angular momentum of the 
reference. Defining the 3-vector n p = 
norm is bounded by 1, this becomes 



L)/(£+|), whose 



i+1 



A+ = / 2 + n p .S 



A; = /2-rvS. 



(4) 
(5) 



Thus, we see that the relational observable is approxi- 
mating a measurement of the spin of the source particle 
along the axis n p = (L)/|(L)|, or rather, of the observ- 
able hp ■ S. In the classical limit, where |(L)| — > £ — > oo, 
the POVM elements become the spectral projectors as- 
sociated to the observable n p • S, and so the measurement 
is perfect. 



B. An operational criterion for the quality of the 
measurement 

Suppose that we set out to measure the the spin of the 
source along the direction n. In other words, we seek to 
measure the projections 



P± = 7 2 /2±n-S. 



(G) 



As we have seen, the measurement {A p } which will ac- 
tually be performed by our method will rather be an 
imperfect simulation of {P^}- There are various opera- 
tional criteria which can be defined to judge the quality 
of this simulation. One such measure [BRST06a | is the 
probability of correctly identifying a maximally polar- 
ized source which is prepared along the directions ±n 
with equal probabilities. This average success probabil- 
ity depends linearly on the inner product n • n p and can 
be written as 



Q 



:./c — 2 (l 



n • n,- 



Other figures of merit could be considered, such as the 
worst-case success probability for identfying the above 
source states, or perhaps another average error in which 
those source states appear with unequal probabilities. It 
is however easy to see that these criteria differ from Q avc 
by 0(l/d) and are thus essentially equivalent. 

There are two distinct properties of the reference that 
contribute to n • n p . On the one hand, the length of n p 
could be less than 1. Indeed, when n and n p are aligned, 
we find tha t n • n n = |nJ. This is the scenario that was 
studied in [BRST06c] | and leads to a "fuzzy" measure- 
ment — a mixing of the + and — measurement outcomes. 
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On the other hand, n and n p may be misaligned, where 
the reference singles out a direction that differs from the 
direction we initially set out to measure. In this case, the 
reference still leads to precise measurements; they just 
don't happen to coincide with what we are interested in. 
An important result that we will establish here is that 
degradation of a reference under repeated uses typically 
lead to misalignm ent rather than fuzzy measurements. 

The authors of BRST06a] defined the longevity of the 
reference to equal the number of times that the reference 
can be used before the above figure of merit drops be- 
low some constant value. They found that the longevity 
of a reference, when used to measure unpolarized source 
particles, scales quadratically with the size of the refer- 
ence. We will see that even if the source particles are 
only slightly polarized, a sufficiently large reference will 
typically have a longevity which scales linearly with its 
size. 



A quantum channel describing the back-action 
on the reference 



In this section we derive exact expressions for the dis- 
turbance experienced by the reference after being used to 
measure a source particle. We will find that just as the 
induced POVM on the source particle depends only on 
the average angular momentum (L) of the reference, this 
back-action can be expressed as a quantum channel, or 
completely positive trace preserving (CPTP) map, which 
depends only on the average angular momentum (S) of 
the source. Since 

J 2 = (i + + §)n+ + (£- \)(t + |)n_, 

a measurement of the total angular momentum is equiv- 
alent to the two-outcome von Neumann measurement 
{n + ,II_}. Suppose the reference and source begin in 
the respective states p and £, after which the above joint 
measurement is performed and the source particle is dis- 
carded. This process causes the reference state to un- 
dergo a discrete time evolution given by the CPTP map 

e c {p) = Tr 8 (n+(p®e)n+ + n_(p®e)n-). 

Using ([3]), we may express this channel as 



i 

2d 2 



+ -Tr s (L-S)(p®C)(L-S) 
+j r ((L.(S»p + p(L.<S»). 



(7) 



It is clear that this expression is independent of the par- 
ticular choice of background coordinates. For a polar- 
ized source, we choose our background coordinate sys- 
tem so that £ commutes with S z , allowing us to write 



£ = 1/2 + 2(S Z )S Z and (S x ) 
Eq. can be written 



(Sy) 



0. In this case, 



1 



1 



2 2d 2 

2(S, 



d 2 



(p 



(L z p + pL z + L + pL- - L-pL + ). 



While not immediately clear from this expression, £^ is 
invariant under rotations about the z axis because the 
only symmetry breaking in JJJ) comes from the polariza- 
tion of the source particles. In the event that the source 
particle is completely unpolarized, this implies that the 
action of the channel £^ is rotationally invariant. In this 
case, £{: reduces to an instanc e of the so-called "spin-i" 
channel introduced in |Rit05a | . A more useful expression 
for can be obtained by straightforward manipulations 
which lead to 



i 

2 

d 2 

1 4 



i - 4(s 2 ; 

2d 2 

L z + (S z ))p(L z + (S z )) 

2{S Z ) 1-2(S Z 
L + pL- + 



(8) 



d 2 



d 2 



-L-pL^ 



When the reference is used to consecutively measure t 
source particles, each of which is in the state £, it will ex- 
perience a disturbance which is described by the channel 
£j = £^ o • • • o £jr where o denotes the usual composition. 
We may think of t as a natural (discrete) time parameter 
for our problem and write the state of the reference after 
it has been used to measure t particles as pit) = £|(jo). 
These dynamics are the central interest of this paper. 

The dynamics of the reference implies a time depen- 
dence of its angular momentum (L(t)) = Tr[p(t)L]. Note 
that because the total angular momentum is a con- 
strained relational observable, the measurement process 
preserves angular momentum. It may however create 
a flow of angular momentum from the reference to the 
source particles. This explains the time dependence of 
(L(£)). This time dependence will in turn induce a time 
dependence of the actual POVM { } which is per- 
formed on the source particle at time step t via Eqs. (|4l5p . 
Thus, the quality of the measurements obtained by re- 
peatedly using the same reference may in general degrade 
over time. 



III. ASYMPTOTIC BEHAVIOR 

We will now analyze the asymptotic state of the refer- 
ence as t ~ * oo. This is an important first step in under- 
standing the general dynamics of a reference. In particu- 
lar, if the asymptotic state of the reference induces good 
measurements of the source, we may conclude that it has 
an infinite longevity. 
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A. Thermodynamical approach 



We may think of the stream of source particles as an 
infinite bath with which the reference interacts. The von 
Neumann entropy of the bath particles is given by 

H(0 = H(i + (S z )) 

where for scalar p, we have written 

H(p) = -plogp- (1-p) log(l -p) 

for the binary Shannon entropy. We may assume that 
the polarization of the source particles is induced by an 
external magnetic field of strength B which is aligned 
with the positive z axis. In this case, the average energy 
of each bath particle is equal to E = —B(S Z ). In units 
where B = 1, we therefore conclude that the bath is at 
inverse temperature 



(3 = 



dH 
~d~E 



log 



1 + 2(S Z ) 
1 - 2{S Z ) 



= 2arctanh(2(S 2 )). (9) 



On physical grounds, we expect the reference to ther- 
malize with the bath after being used to measure a large 
number of particles. The asymptotic partition function 
is thus anticipated to be 



Z = Trexp(-/3L Z ) 



sinh(/3d/2) 
sinh(/3/2) ' 



(10) 



leading to the asymptotic reference state = 
4 exp(— (3L Z ). In the next subsection, we prove that this 
corresponds to the asymptotic state by showing that it 
is the unique fixed point of the quantum channel 

Meanwhile, let us investigate the usefulness of this 
asymptotic state to serve as a reference. Clearly, the an- 
gular momentum (L) of the reference in its asymptotic 
state is along the z direction. Thus, if the original pur- 
pose of the reference was to measure the spin of the source 
particles along an axis that differs significantly from z, we 
conclude that it will fail to do so asymptotically simply 
because it will be misaligned. 

On the other hand, if the reference is meant to measure 
the spin of the source particles along their axis of polar- 
ization, then it may continue to serve this purpose in its 
asymptotic state. From our discussion in Sec. Ill Bl we 
know that the figure of merit is in this case given by the 
length of the vector n p = (L)/(£+ ~). The expectation of 
the z component (L z ) of the asymptotic reference's angu- 
lar momentum, as a function of the inverse temperature 
(3, is given by 



1 dZ 
Z~d(3 

sinh(/3/2) d sinh(/3d/2) 
sinh(/M/2) d(3 sinh(/3/2) 
(d coth(/3d/2) - coth(/5/2)) /2. 
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FIG. 1: Normalized polarization of the reference {L z )/£ ther- 
malized with a source of polarization (S z ) for the various ref- 
erence sizes I — 20, 40, and 80. The transition sharpens as 
the reference's size increases, with a width of roughly l/£, and 
to first order, we find that the slope of (L z ) is either £ or l/l 
accordingly. In the macroscopic limit, any source polarization 
induces a full reference polarization. 



Using ©, we may rewrite this expression in terms of the 
source polarization (S z ) to obtain 



(L z ) = -coth (darctanh(2(S z ))) 



HS Z 



Note that this expression is an antisymmetric function 
of (S z ). In particular, the asymptotic reference state is 
well-defined when (S z ) = by continuity. In Fig. [T] we 
have plotted the normalized quantity (L z }/£ for various 
reference sizes £. From this, we can compute the figure 
of merit |nJ. Observe that the case of an unpolarized 
source (i.e. (S z ) w 0) is quite singular and leads to a 
completely different physical behavior of the reference 
than in the generic case. For the case of a vanishing 
source polarization in which (S z ) <C l/l, the figure of 
merit vanishes: 



Ins, 



O 



(Sz 

I 



On the other hand, for a source with a finite amount 
of polarization, where (S z ) > £~ a with < a < 1, we 
instead find that the figure of merit approaches 1: 



n. 



> 1 - 



1 



4 fi 



B. Proof of thermalization of the reference 

To check that the reference asymptotically reaches the 
thermal state discussed in the previous section, one can 
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resort to a brute- force calculation that exp(— 0L Z ) is a 
fixed point of Another instructive way to see this 
begins by writing p in the basis \m) which diagonalizes 
L z . By inspecting Eq. [5J we see that £^ does not mix the 
various diagonals of p, i.e. 



Ss(\m)(m'\) 



1 

E 

a=-l 



P(m+a, m '+a)\( m ,m>)\m+a)(m' + a\ 



for nonnegative numbers P( m+a ,m'+a)\{m,m') satisfying 
1 

^ -^(m+a,m'-j-a) | (m.m') — 1- 
<j=-l 

One may further check that this bound is saturated 
only when m = m' . Therefore, when m ^ m! , these 
weights correspond to subnormalized conditional proba- 
bilities which asymptotically annihilate the off-diagonal 
terms of p. Abbreviating 

Pm+a^n — P(m+a : m+a) | (m,m) ) 

we may express the action of £^ on diagonal states (those 
which commute with L z ) as 

l 

£ £ (\m){rn\) = ^2 P m+a \ m \m + a){m + a\. 

The conditional probabilities can be directly evaluated 
using Eq. [51 obtaining 

1 ± 2(S Z 



P 



mzbllm 




and P, 



1 — Pm+i\m — Pm-i\m- This defines a Markov 



chain whose stationary state = 5Zm ^"l m )( m l ^ s 
found by demanding the detailed balance equation 

PmPrn+l\m — Pm+lPm\m+l ) 

which yields the unique solution 

Pm(X i}rW)) =exp ^- 

This proves that the asymptotic quantum partition func- 
tion of the reference is indeed Z = Trexp(— (3L Z ). 

IV. THERM ALIZATION TIME 

So far, we have derived expressions for the asymptotic 
state of the reference and have characterized its ability 
to measure the spin of the source particles in this state. 
We will now analyze how much time it takes for the ref- 
erence to reach this asymptotic state. In circumstances 
where the asymptotic state is not a useful resource - 
i.e. when it becomes misaligned or fuzzy — this ther- 
malization time will set a fundamental upper bound on 
the longevity of the reference. 



A. A lower bound on the thermalization time via 
conservation laws 

Here, we give a straightforward derivation of a lower 
bound on the thermalization time, using only basic sym- 
metry arguments. Since measuring the total angular mo- 
mentum of the combined system is a rotationally invari- 
ant process, it preserves angular momentum. To signif- 
icantly disturb the state of a reference with large an- 
gular momentum requires "kicking it" with comparable 
amount of angular momentum. When each source parti- 
cle is in a state with angular momentum (S z ), the total 
angular momentum of the source particles is a binomial 
random variable. Therefore, the typical states of t source 
particles will have a total angular momentum given by 



5> 



typical 



= t(s,)±Ji-(s,)*tf 



In the case of an unpolarized source (((Sz)) < 1/^), 
the physical requirement to disturb the state of the ref- 
erence leads to a number of measured particles t > £ 2 ; 
the disturbance is due to fluctuations of the source's 
angular m o mentum. This explains the main results of 
[BRST06at lBRST06b| . On the other hand, when the 
source is sufficiently polarized (|(<S Z )| > l/£), the distur- 
bance is no longer dominated by fluctuations but rather 
by the bias, or drift, leading to a longevity t > £/2(S z ). 

This argument is not restricted to the particular SU (2) 
model we have been studying here; it should extend to 
arbitrary reference systems as well. In particular, when 
the source particles are in random states, the expecta- 
tion value of any traceless local operator will grow as 
the square-root of the number of particles. This yields a 
strict lower bound on the longevity of the reference. 



B. An upper bound on the thermalization time 



We have demonstrated in Sec. lIII El how the map £j re- 
duces to a Markov chain on density operators diagonal in 
the L z basis. Thus, we can use standard tools to analyze 
the expected hitting time for any state. The expected 
hitting time of a state \£, m) is the average time it takes 
the system to reach this state. 

Consider the case where the reference begins in a co- 
herent state which points in the — z direction, i.e. it is 
initially in the pure state \£, —£), and the source is maxi- 
mally polarized, i.e. (S z ) — ^. Intuitively, we expect this 
"antiparallel" setting to be a worst case scenario, and so 
to yield a fundamental lower bound on the longevity of 
the reference. 

Substituting the value of (S z ) into Eq. (fTT]l . we see 
that at every time step, the angular momentum of the 
reference will either increase by ^ , or will otherwise stay 
the same. In this case, the unique steady state is \£, £) as 
expected. The average time T required to arrive at this 
state starting from the state \£, —£) can be computed by 
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adding the average time required to go from state \£, m) 
to \£,m+l) for m = -£,-£ + !,...,£- 1: 



^--^ -i _ (2m+l\ 2 
,l- £ x 

~ d / t - 

= d\og{d- 1) 

where e = 2/d. If the reference is in the fixed point, 
it gives the incorrect result with certainty: this is not 
surprising since it started pointing in the — z direction 
and ended up pointing in the +z direction. Thus, T 
gives an upper bound on the longevity of the reference. 

Note that the logarithmic factor is due to our con- 
vergence criterion becoming more strenuous as we in- 
crease £. Had we instead computed the average time to 
reach a state with (L z )/£ > 1 — e starting in a state 
with (L z ) 1 1 < 1 — e where e is a constant indepen- 
dent of d, we would have obtained a linear hitting time 
T = dlog(- — 1). Since this matches the lower bound 
obtained in the previous section, we conclude that this 
scaling is asymptotically tight. Hence, we conclude that 
the longevity of the reference in this worst case is upper 
bounded by a linear function of the size of the reference. 

For more general settings, this thermalizing time can 
be evaluated using the spectral gap theorem. Given the 
stochastic Markov matrix P mm > defined by Eq. (fTT|) . we 
construct the matrix W = p(oo)~ 1 I 2 P p(oo) 1 I 2 . The con- 
vergence time is equal to the inverse of the gap between 
the two largest eigenvalues of W. We have found numer- 
ically that this gap is given by A w 2(S z )/£ as pictured 
in Fig. [21 Thus, the time required to reach the fixed 
state is t oc £/2(S z ), once again in agreement with the 
predictions of the previous section. 




100 „ 200 300 



FIG. 2: Gap inverse (A^ 1 ) of the matrix W as a function 
of the reference size £ and for different values of (S z ) = 
0.05,0.1,0.2, and 0.3. Dashed lines indicates the estimate 



V. NON-EQUILIBRIUM STATE 

We have studied the asymptotic state of the reference 
and characterized the time it takes to get there. Our fi- 
nal task is to study the thermalization process itself — 
to describe the dynamics of the reference as it reaches 
its equilibrium state. More specifically, since the expec- 
tation (L}/£ fully characterizes the operational behavior 
of the reference, our aim is to derive an equation for the 
time evolution of this quantity. As we did before, we 
choose our coordinates so that the source's polarization 
is along the z axis and the initial orientation of the ref- 
erence is at an angle 9 from the z axis in the z — x plane, 
or in other words (L y ) = 0. By symmetry, we know that 
(L y ) will remain for all times. Thus, we can parame- 
terize n by its length r and its angle 9 away from the z 
axis. Figure [3] displays the thermalization trajectories of 
the reference for different reference sizes I and different 
initial coherent states |#o) = exp{— i9oL y }\£, I). 




(L x )/£ (L x )/£ (L x )/£ 

FIG. 3: Phase space trajectories ({L x (t)}/£, (L z (i))/£) of the 
equilibration process for different initial states \6o), and dif- 
ferent reference sizes £ — 20, 40, and 80. The dashed 
line represents the manifold of coherent states, for which 
(L x ) 2 + (L z ) 2 = I 2 . The source is maximally polarized 
(S z ) = 1/2. The loss of polarization of the reference dur- 
ing equilibration is suppressed as its size increases. In the 
macroscopic limit, the reference remains in a coherent state 
throughout the process, except for the singular initial state 
0o = 7T. 



To derive an equation of motion, it is convenient to 
express the CPTP map £(. from Eq. ([7]) acting on the 
reference using a rotated basis of operators L , given by 

(L 9 X , L e y , L e z ) = (cos 6L X - sin 9L Z , L y , sin 9L X + cos 9L Z ) . 

In particular, we choose 9 so that (L x ) = (Ly) = and 
(L e z ) = r£ for some < r < 1. With some straightforward 
manipulations, we can express this map in the Kraus 
form £^(p) = Y] a E a pE\ where the four Kraus operators 
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E a are given by 

^{v^ 2 + l-4(S 2 ) 2 , i2y/l-l(S z )*L e y 
2L e z + iA{S z ) sin 9L e y + 2(S Z ) cos 6, 
2L X + iA(S z ) cos 9L e y - 2{S Z ) sin o\. 

To make further progress, observe that when acting on 
the state of the reference, L e z ~ £r and L e x ,L e y ~ V?- 
Expanding in powers of 1 and keeping only first-order 
terms, we obtain the approximation 



r(S z ) 

£ t {p) k p + i—^- sin 0[L y ,p\. 



(12) 



Observe that the second term induces a rotation by an 
angle rzsm9/£ about the y axis. This approximation 
leads to the following equation for the angle 9{t) relative 
to the z axis 



smf 



o(i/e 2 ). 



(13) 



By numerical analysis which is summarized in Fig. [3J 
we see that in the macroscopic limit, the polarization of 
the reference remains close to 1 throughout the entire 
trajectory, provided that the initial polarization of the 
reference is sufficiently far from being anti-parallel with 
that of the source. Assuming that r(t) = 1 throughout 
the entire trajectory, we may thus solve Eq. (|13[) to obtain 
a semi-classical trajectory: 

e sc (t) = 2arccot(cot(6»o/2)e <s ^* /£ ). 

This solution is compared with the exact numerical eval- 
uation for various reference sizes in Fig. 0J We note that 
this der ivation jus tifies the ad hoc semi-classical model 
used in BRST06b] . 



VI. DISCUSSION 

We have analyzed the behavior of a quantum gyroscope 
which is modeled as a large spin and used to measure se- 
quences of spin-i source particles from a fixed ensemble. 
We showed that the measurement which is actually per- 
formed on the source depends only on the angular mo- 
mentum of the reference. We then derived expressions 
for the quantum channel which describes the back-action 
on the reference and showed that its unique fixed point 
is a thermal state. It is also possible to characterize the 
asymptotic state of the reference as 



p(oo) ex n sym ^n ; 



sym ; 



where n sym is the projector onto the (2£+ 1 dimensional) 
symmetric subspace of (C 2 )® . This expression is in- 
sightful when the reference is a composite system made 
up of 2£ spin-i particles, as one can imagine that the 




FIG. 4: Orientation of the polarization of the reference rela- 
tive to the z axis as a function of time. The source's polariza- 
tion is non maximal (S z ) = 1/4 and the reference is initially 
in the coherent state 8o = 157r/16. The dotted line is the 
prediction obtained from Eq. (|13p . Solid lines represent the 
exact numerical solutions for I = 20, 40, and 80. As the size 
of the reference increases, the exact trajectories approach the 
semi-classical solution 9sc(t). 



particles of the reference get "exchanged" with those of 
the source, while remaining in the spin-£ representation. 

We then estimated the amount of time it takes to ap- 
proach the thermal state. For a polarized source, we saw 
that the thermalization time scales linearly with the size 
of the reference. The implications of this for the longevity 
of the reference are two-fold. If the source and reference 
begin in alignment, the longevity will be infinite for a 
sufficiently large (but still finite) reference size. This is 
because the asymptotic state will still have most of its 
angular momentum pointing in the original direction. In 
all other cases with a polarized source, the reference will 
reorient itself to become aligned with the source after a 
time that scales linearly with its size. 

We also g ave a physic al mechanism describing the main 
result from |BRST06al |: for an unpolarized source, the 
thermalization time scales quadratically with the refer- 
ence size because the dynamics is dominated by fluctu- 
ations rather than by drift. Therefore, its longevity en- 
joys the same scaling. However, we saw that this type 
of behavior is quite non-generic; unless the angular mo- 
mentum of the source is exactly zero, a sufficiently large 
(but finite) reference will asymptotically retain a near 
maximal amount of angular momentum. The quadratic 
scaling, on the other hand, requires a precisely tuned 
source, reminiscent of the effort needed to ensure that 
a statistical system remains in the vicinity of its critical 
point. 

Finally, we initiated the analysis of the non-equilibrium 
dynamics of the reference, giving a semiclassical equation 
for the time-evolution of the reference. In particular, we 
found that if the reference begins in a coherent state, 
it will essentially remain in a coherent state throughout 
its evolution, provided that the source particles are not 
anti-parallel to the initial state of the reference. Here, 
the reference simply rotates itself toward the final state 
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in an essentially deterministic fashion. Thus, while it 
loses its ability to perform the measurement for which 
is was originally intended, it will always be reasonably 
good at performing some particular measurement (which 
is predicted by our equation) at any given time. For an 
even better understanding of the semi-classical dynam- 
ics, it would be beneficial to also have an equation for 
the time evolution of the length r of the reference, or at 
least to obtain an analytic bound on its minimal value 
throughout its evolution. 

As mentioned in the introduction, our analysis is rel- 
evant for the realization of quantum information pro- 
cessors. For micro-fabricated architectures, some com- 
ponents of the measurement device can be subject to 
quantum fluctuations and will thus feature effects sim- 
ilar to those analyzed in this paper. Most sc hemes fo r 
fault-tolerant quantum computation (see e.g. Sho96a| ) 
utilize repeated measurements for syndrome extraction 
and post-selected state distillation. Moreover, the mea- 
sured particles will typically exhibit a strong bias (e.g. 
syndromes are more likely to be than 1), inducing a 
rapid degradation of the reference. On the other hand, 
it may also be possible to utilize such a source bias for 
stabilizing a reference, allowing sufficiently accurate mea- 
surements to be performed throughout an arbitrarily long 



computation. 

One could imagine extending the analysis put forth in 
this paper in a number of ways. For example, the refer- 
ence could be used to measure the angular momentum of 
particles in higher representations of SU(2). Addition- 
ally, more complicated Lie groups could be substituted; 
repeating our analysis in such cases would inv olve even 
further generalizations of channels studied in [Rit05a|. 
As many of our arguments are based on general physi- 
cal principles, we expect that the same scaling properties 
derived here would apply in those cases as well. 
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